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GPB SIMULATION DEVELOPMENT FINAL REPORT 


1. Introduction 

This report documents the simulation of the Lockheed Martin designed Gravity 
Probe B (GPB) spacecraft developed tool by bd Systems Inc using the TREETOPS 
simulation. This study quantifies the effects of flexibility and liquid helium slosh on GPB 
spacecraft control performance. The TREETOPS simulation tool permits the simulation 
of flexible structures given that a flexible body model of the structure is available. For 
purposes of this study, a flexible model of the GPB spacecraft was obtained from 
Lockheed Martin. To model the liquid helium slosh effects, computational fluid 
dynamics (CFD) results 1 were obtained from Dr. R. J. Hung at the University of 
Alabama in Huntsville (UAH), and used to develop a dynamic model of the slosh 
effects. The flexible body and slosh effects were incorporated separately into the 
TREETOPS simulation, which places the vehicle in a 650 km circular polar orbit and 
subjects the spacecraft to realistic environmental disturbances and sensor error 
quantities. In all of the analysis conducted in this study the spacecraft is pointed at an 
inertially fixed guide star (GS) and is rotating at a constant rate about this line of sight. 


2. GPB Simulation Components 

The GPB TREETOPS simulation is made up of the spacecraft control sensors and 
actuators, the attitude and translation controllers, environmental disturbances, a slosh 
model of the liquid Helium, a flexible body spacecraft dynamics model, and a science 
gyro model. Each of these models is discussed in the following subsections. 


2.1. Spacecraft Configuration 


2.1.1. Sensors 

The spacecraft sensors modeled in this simulation are the science telescope, drag 
free sensor, control gyros, and the roll star trackers. The sensor errors specified in this 
section were taken from Lockheed Martin 2 . 

The science telescope is modeled using the TREETOPS star tracker sensor. This 
sensor provides a line of sight error about the spacecraft pitch and yaw axes during 
guide star valid (GSV) phases. The science telescope model for each axis is shown in 
Figure 2.1. 1-1. Science telescope error is modeled by a zero mean gaussian noise 
with a standard deviation of 22.36 milliarcseconds that is added to the true pointing 
error. The resulting signals are then quantized using a least significant bit weight of 
0.25 milliarcseconds. 


1 
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Figure 2.1 .1-1 Science Telescope Model 

Noise 



The control gyros are modeled using the TREETOPS rate gyro sensor. These 
sensors provide the angular rate with respect to inertial space about each of the three 
body axes. The control gyro model for each axis is shown in Figure 2.1. 1-2. For the 
pitch and yaw axes, a gaussian noise with a mean drift of 0.003 arcsec/second and a 
standard deviation of 0.002236 arcsec/second is added to the true rate, then the 
signals are quantized using a least significant bit weight of 1 .3 milliarcseconds. For the 
roll axis, a gaussian noise with a mean drift of 0.003 arcsec/second and a standard 
deviation of 0.002236 arcsec/second is added to the true rate, after which the signals 
are quantized using a least significant bit weight of 37.5 milliarcseconds in low range 
mode or 375 milliarcseconds in high range mode. During the guide star invalid (GSI) 
phase, the control gyro derived rates are propagated using quaternions to produce the 
pitch and yaw line of sight errors. At the end of each guide star valid (GSV) phase, the 
science telescope signal is used to update the initial attitude to be used in the attitude 
propagation during the guide star invalid phase. 

Figure 2.1.1-2 Control Gyro Model 


Noise 



The roll star tracker is modeled using the TREETOPS integrating rate gyro sensor. 
This provides the roll angle. The roll star tracker model is shown in Figure 2.1. 1-3. A 
zero mean gaussian noise with a standard deviation of 5.0 arcseconds is added to the 
true roll angle. The resulting signal is then quantized using a least significant bit weight 
of 0.5 arcseconds. The roll error is computed by subtracting this signal from the 
commanded roll angle. 


2 
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Figure 2.1. 1-3 Roll Star Tracker Model 


Noise 



The drag free sensor is modeled using the TREETOPS position vector sensor. 
This provides a relative position vector from the spacecraft to the drag free proof mass 
in the inertial frame. The inertial frame relative position vector is converted to the body 
frame. For each translational body axis, a zero mean gaussian noise with a standard 
deviation of 3.354 nanometers is added to the true displacement, then the signals are 
quantized using a least significant bit weight of 1 nanometer. 

Figure 2.1 .1-4 Drag Free Sensor Model 


Noise 



2.1 .2. Actuators 

The 16 thrusters are modeled using 16 of the TREETOPS reaction jet actuators. 
Each thruster is modeled as shown in Figure 2. 1.2-1. The commanded thrust for each 
thruster is quantized using a least significant bit weight of 0.0025 milliNewtons. Then 
each thrust signal is delayed one control cycle (0.1 seconds) and a hysteresis of 0.005 
milliNewtons is applied 3 . Finally, a separate zero mean gaussian noise with a standard 
deviation of 0.0559 milliNewtons is added to each thruster. 

Figure 2.1 .2-1 Thruster Model 


Noise 



Quantization Delay Hysteresis/Umiter 
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2.2. Environmental Disturbances 

The environmental disturbances implemented in the GPB TREETOPS simulation 
include the cryoperm shield magnetic torque, the aerodynamic force and torque, and 
the gravity gradient force and torque. 

2.2.1 . Magnetic Disturbance 

The cryoperm shield magnetic moment is calculated as 

[176320 0 0 

M M( = 0 160320 0 B b 

0 0 2160000 

where B b is the Earth’s magnetic field vector in the spacecraft body frame. The 
cryoperm shield magnetic torque vector is calculated as 

7 ctjo ~ x B b • 

The cryoperm shield magnetic torque for one orbit is shown in Figure 2.2-1 and is 
applied at the spacecraft center of mass. 

2.2.2. Aerodynamic Disturbance 

The aerodynamic force and torque vectors in the spacecraft body frame are calculated 
as follows 

0.0025 sin((0 (> r +7c/2)sin(co r f - it/2) 

= 0.0025 sin(o) a t + 7t/2)sin(a) r t) 

0.0035 sm(co 0 r+rc) 


0.002 sin((o e r)sin(2co o r)sin(a) r t + jc) 

= 0.002 sin((o < ,r)sin(2o) 0 r)sin(a),r + Ji/2) 

0.0024 sin (co 0 t + 7t/2)sin(2co o r) {1 + 0.06 sin(4cfl ,f )} 

where co, is the orbit rate and co, is the roll rate. The aerodynamic force and torque for 
one orbit are shown in Figures 2.2-2 and 2.2-3, respectively, and are applied at the 
spacecraft center of mass. 


2.2.3. Gravity Gradient Disturbance 
The gravity gradient force and torque vectors are computed as 

R 0 xIR 0 + 

, 3|X . >2 
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where p. is the gravitational constant, m K is the spacecraft mass, is the vector from 
the center of mass to the proof mass, R 0 is the vector from the center of the Earth to the 
spacecraft center of mass and I is the spacecraft inertia dyadic, N is the north vector, 
and h is the gravitational coefficient 0.001 08263. The gravity gradient force and torque 
vectors are shown in Figures 2.2-4 and 2.2-5, respectively, and are applied at the 
spacecraft center of mass. 
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Figure 2.2-1 Cryoperm Shield Magnetic Disturbance Torque 


TREETOPS GPB SIM: Magnetic Toque (roll rate = 0.1 rpm) 
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Figure 2.2-2 Aerodynamic Disturbance Force 


TREETOPS GPB SIM: Aerodynamic Force (roll rate = 0.1 rpm) 
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Figure 2.2-4 Gravity Gradient Disturbance Force 


TREETOPS GPB SIM: Gravity Gradient Force (roll rate = 0.1 rpm) 
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2.3. Control System Operation 


2.3.1. Translational Control 

The GPB translational controller* used in the TREETOPS simulation is shown in 
Figure 2.3.1-1. This same controller is used for all three axes. The translational 
position error (i.e. the distance between the spacecraft and proof mass along each 
body axis) comes from the drag free sensor. This error is converted from the body 
frame to the nadir frame for control computations. A PID controller is used with a limiter 
on the integral loop to produce a force command. This command is filtered by a 
second order filter which has a peak at the roll frequency to attenuate disturbances at 
the roll frequency. The output of the filter is transformed back to the body frame and 
then limited before being sent to the thruster selection logic. The translational controller 
gains and limits are shown in Table 2.3.1 *1. 

Table 2.3.1-1 Translational Controller Gains and Limits 


LIMITS 


1 (micrometers) 

0.5 

2 (milliNewtons ) 

5.0 



GAINS 


K p (1 / sec 2 ) 

2.44 

K r ( sec ) 

1.2435 

K ( 1 / sec ) 

0.605 



Bandwidth ( Hz ) 

0.5 
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FIGURE 2.3.M Translational Controller 
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2.3.2. Attitude Control 

The GPB attitude controller 4 used in the TREETOPS simulation is shown in Figure 
2.3.2-1 . This same controller is used for all three axes. The attitude error for pitch and 
yaw comes from the science telescope during Guide Star Valid phases and from the 
propagated control gyro rates during Guide Star Invalid phases. The attitude error for 
roll comes from the roll star trackers. The attitude rate errors for the three axes come 
from the control gyros. A PID controller is used with a limiter on the proportional, 
integral, and derivative loops, as well as on the filtered sum. The control command is 
filtered by a second order filter which has a peak at the roll frequency to attenuate 
disturbances at the roll frequency. The output of the attitude controller is a torque 
command for each axis which is sent to the thruster selection logic. The attitude 
controller gains and limits are shown in Table 2.3.2-1 . 

Table 2. 3.2-1 Attitude Controller Gains and Limits 


LIMITS 

Pitch/Yaw 

Roll 

1 (arcsec ) 

1.0 

60.0 

4 (arcsec ) 

0.25 

20.0 

5 (arcsec ) 

0.05 

4.0 

6 (arcsec ) 

0.25 

20.0 

7 (arcsec ) 

0.25 

20.0 




GAINS 



Kp ( 1 / sec 2 ) 

0.482 

0.054 

K r ( sec ) 

1.98 

7.26 

K, ( 1 / sec ) 

0.25 

0.0023 




Bandwidth ( Hz ) 

0.15 

0.06 
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FIGURE 23.2-1 Attitude Controller 
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2.3.3. Thruster Selection Logic 

The thruster selection logic determines the thrust command to each of the sixteen 
thrusters based on the commanded force vector from the translation controller, the 
commanded torque vector from the attitude controller, and the helium flow command. 
The flow diagram of the thruster selection logic is shown in Figure 2.3.3-1 . 


Figure 2 .3.3-1 Thruster Selection Logic 

Compute the 8x6 pseudo-inverse of the 6x8 transformation matrix. A, 
of thruster pair forces (8) to control forces and torques (6) 

A' 1 = A T (A A 1 )* 1 


Compute the thruster pair forces (8) from the control forces and 
torques using the pseudo-inverse of A 


Fpa® — 



Tq - FpaIR i 


II 

o 

b 

+ 

II 

o 

o 


lq+i - f pAiRi 


If the helium flow command is greater than the helium flow needed for 
control, add the residual mass flow to each thruster command, T G to 
satisfy the flow command requirements 

Tc,=T cl +*fm,I,Ml- XtJ 

V i=ltol6 ) 
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2.4. GPB Helium Slosh Model 

The slosh of the Helium in the dewar of the GPB spacecraft is modeled in the 
TREETOPS simulation as shown in Figure 2.4-1. Four slosh masses are placed 
symmetrically about the dewar center of mass, two on the x axis and two on the y axis. 
Each of these masses is free to move in all three coordinate directions. The distance 
from the vehicle centerline (z axis) to the connect point of each mass is determined by 
the spin rate of the spacecraft and the stiffness of the spring connecting each slosh 
mass to the main body. The initial displacement of the springs between the main body 
and slosh masses are set to start the simulation in an equilibrium condition. Table 2.4- 
1 shows the slosh input data that was chosen to match the results of Dr. Hung 1 . 

Table 2.4-1 Helium Slosh Model Input Data 


Mass of each slosh body 

71 .9 kq 

Number of slosh masses 

4 

Dewar tank fill level 

80% 

Distance from centerline to slosh mass 

0.35 m 

Initial sprinq displacement 

0.22778264 m 

Sprinq Stiffness (x, y, z) 

0.02, 0.02, 0.04 N/m 

Sprinq Dampinq (x, y, z) 

0.5, 0.5, 0.01 N/(m/s) 

Spacecraft spin rate 

0.1 rpm 

Distance from spacecraft cm to slosh cm (z) 

-0.1 m 
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2.4.1. GPB Helium Slosh Model Comparison Results 
The slosh model results from the simulation were compared to the results from Dr. 
Hung of the University of Alabama at Huntsville (UAH)'. Figure 2.4.1-1 shows the 
comparison of the helium fluid center of mass motion due to the aerodynamic 
force/torque disturbances for one orbit. Figure 2.4.1-2 shows the comparison of the 
helium fluid center of mass motion due to the gravity gradient force/torque disturbances 
for one orbit. Figure 2.4.1 -3 shows the comparison of the helium fluid center of mass 
motion due to the cryoperm shield magnetic torque disturbance for one orbit. The GPB 
TREETOPS simulation (bd Systems) results of fluid center of mass motion compare 
well with Hung’s results for the aerodynamic and gravity gradient disturbances, with the 
only significant difference being a different high frequency component in the z axis (roll) 
results. The comparison of fluid center of mass motion due to cryoperm shield 
magnetic torque was not as close as with the other disturbances. Retuning of the 
model was not accomplished, however, because it was deemed most important to 
match the fluid motion caused by the largest disturbances i.e. aerodynamic and gravity 
gradient effects. 
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Figure 2.4.1-1 Helium Fluid Center of Mass Motion Due To The Aerodynamic 

Force/Torque Disturbances 



TJ 
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Figure 2.4.1-2 Helium Fluid Center of Mass Motion Due To The Gravity Gradient 

Force/Torque Disturbances 
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Figure 2.4.1-3 Helium Fluid Center of Mass Motion Due To The Cryoperm Shield Magnetic 

Torque Disturbance 
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2.5. GPB Flexible Body Model 

The TREETOPS GPB flexible body spacecraft model was generated from an 
on-orbit structural dynamics model* supplied by Lockheed Martin. A graphical 
representation of this model is shown in Figure 2.5-1. This model was supplied to 
MSFC in an l-DEAS Master Series 2.1 format It was converted by MSFC to a MSC 
Nastran universal input format and given to bd Systems for processing. The model 
represents the beginning of mission deployed configuration. The mass properties are 
shown in Table 2.5-1 and the node locations used in the TREETOPS model are shown 
in Table 2.5-2. The closest node in the model to any of the science gyro locations was 
at grid 10950, whose location was between science gyros 2 and 3. The modal 
frequencies are shown in Table 2.5-3. These first eight flexible body modes were used 
in the TREETOPS simulation and are all the flexible body modes below 13 Hz. 

Figure 2.5-1 GPB Flexible Body Spacecraft FEM Drawing 
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Table 2.5-1 GPB Finite Element Model (FEM) Mass Properties 



FEM 

Mass (kg) 

3190 

Ixx (kg-m 2 ) 

4320 

lyy (kg-m 2 ) 

4418 

Izz (kg-m 2 ) 

2760 

Cgx (m) 

0 

Cgy (m) 

0 

cg&M 

0 


Table 2.5-2 GPB Flexible Body Model Node Locations 


Node 

Node No. 


Y(m) 

_z(m) 

Center of Mass 

13001 

0.0 

0.0 

0.0 

Control Gyro 

40002 

0.19050 

1.08457 

0.60196 

Proof Mass 

10949 

0.0 

0.0 

-0.04810 

Thruster Pod 1 

921230 

-1.28651 

0.0 

2.547830 

Thruster Pod 2 

921231 

1.28651 

0.0 

2.547830 

Thruster Pod 3 

921229 

-1.28651 

0.0 

-1.86217 

Thruster Pod 4 

921228 

1.28651 

0.0 

-1.86217 


10950 

0.0 

0.0 

-0.37873 


Table 2.5-3 GPB Flexible Body Model Modal Frequencies 


Mode No. 

Frequency (Hz) 

1-6 

0.0 

7 

1.32 

8 

1.63 

9 

1.93 

10 

2.09 

11 

2.28 

12 

2.72 

13 

3.78 

14 

4.57 
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2.6. GPB Science Gyro Model 

A simple science gyro model was developed to compute the north-south and east- 
west Newtonian drifts in the science gyro. A graphical depiction of the science gyro 
model is shown in Figure 2.6-1. According to Kasdin and Gauthier*, the Newtonian 
drifts are caused by two sources, electrostatic support dependent torques and 
electrostatic support independent torques. The dependent torques include the earth's 
gravity gradient force acting on rotor mass unbalance and on rotor oblateness (not 
modeled). The independent torque is that produced by earth's gravity gradient acting 
directly on the oblate science gyro rotor. The science gyro center of gravity is assumed 
to be displaced from its center of geometry by 2 pin. The electrostatic forces are 
applied at the science gyro center of geometry. Only forces perpendicular to the 
angular momentum vector are assumed to produce noncyclic gyro drift 

2.6.1. GPB Science Gyro Comparison Results 

The science gyro drift angles due to spherical earth gravity gradient force acting 
through the rotor mass unbalance are shown in Figure 2.6.1-1, with a negligible east- 
west drift of -2.6x1 O'* macs/year. The science gyro drift angles due to an oblate earth 
gravity gradient force acting through the rotor mass unbalance are shown in Figure 
2.6.1 -2, with an east-west drift of -1 .77 macs/year. The science gyro drift angles due to 
an earth gravity gradient torque, including J2 effects, on an oblate rotor are shown in 
Figure 2.6.1 -3, with an east-west drift of 0.013 macs/year. 
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FIGURE 2.6-1 Science Gyro Model 
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Figure 2.6.1-1 Science Gyro Drift Angles Due To Spherical Earth Gravity Gradient Force 

Acting On Rotor Mass Unbalance 


SG NORTH-SOUTH DRIFT ANGLE (Fgg) 



x in'* SG Gravity Gradient Force 
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Figure 2.6.1 -2 Science Gyro Drift Angles Due To An Oblate Earth Gravity Gradient Force 

Acting On Rotor Mass Unbalance 
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Figure 2.6.1 -3 Science Gyro Drift Angles Due To An Earth Gravity Gradient Torque, 

Including J2 Effects, On An Oblate Rotor 
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3. GPB Simulation Performance Results 

The purpose of the performance evaluation was to determine the spacecraft 
pointing control response, in both the guide star valid and invalid phases, the roll 
control response, and the translation control response. All significant disturbances 
affecting the vehicle were evaluated, including aerodynamic disturbances, gravity 
gradient disturbances, and cryoperm shield magnetic disturbances. All sensor and 
actuator error characteristics were included in the performance evaluation. Individual 
runs were made to evaluate the spacecraft’s ability to meet the budgeted control 
requirements. A summary of the performance runs is shown in Table 3-1. The 
performance responses for the flexible body runs accomplished at .3 rpm and .5 rpm 
are shown in Figures 3-1 through 3-38. 


Table 3-1 Performance Run Summary 


RIGID BODY @ 0.3 rpm ROLL RATE 

GSV LOS 
Error 

(marcs) 

■ 

■ 

I 

HHBH 

ENVIR (AERO+GRAV GRAD+CRYO) 

0.33 

2.18 

0.96 

0.36 

THRUSTER ERRORS 

9.57 

8.63 

0.018 

7.41 

ROLL STAR TRACKER 
NOISE/QUANTIZATION 

0.55E-1 1 

1.0E-11 

0.028E-6 

1.0E-10 

SCIENCE TELESCOPE 
NOISE/QUANTIZATION 

7.82 

21.37 

0.54 

0.31 

CONTROL GYRO DRIFT 

0.0051 

191.9 

0.047E-3 

0.0033 

CONTROL GYRO NOISE 

1.29 

32.93 

0.0015 

0.079 

CONTROL GYRO QUANTIZATION 

0.17E-10 

0.22E-1 0 

0.002 

0.45E-1 1 

DRAG FREE SENSOR NOISE 

0.15E-4 

0.16E-4 

0.28E-6 

3.10 

DRAG FREE SENSOR QUANTIZATION 

0.50E-1 1 

0.60E-1 1 

0.30E-6 

0.13E-1 1 


FLEXIBLE BODY @ 0.3 rpm ROLL RATE 

GSV LOS 
Error 

(marcs) 



1 

1H 1 

ENVIR. (AERO+GRAV GRAD+CRYO) 

0.33 

3.66 

0.96 

0.36 

THRUSTER ERRORS 1 

9.64 

8.65 

0.018 

7.37 

ROLL STAR TRACKER 
NOISE/QUANTIZATION 

0.28 

■B 

0.54 

■fil 

SCIENCE TELESCOPE 
NOISE/QUANTIZATION 

7.90 

22.88 

0.13E-4 

0.32 

CONTROL GYRO DRIFT 

0.0054 

192.9 

0.47E-3 

! 0.0034 

CONTROL GYRO NOISE 

1.31 

33.05 

0.0015 

0.082 

CONTROL GYRO QUANTIZATION 

0.77 

1.96 

0.002 

0.61 

DRAG FREE SENSOR NOISE 

0.03 

1.76 

0.34E-4 

3.10 

DRAG FREE SENSOR QUANTIZATION 

0.0016 

1.77 

0.20E-5 

1.23 
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Summary Runs 





CONTROL (ENVIRONMENT & THRUSTERS) 

9.52 

8.97 

0.96 

7.24 

DETERMINATION (ALL SENSORS) 

7.90 

215.66 

0.70 

3.44 

TOTAL 

(ENVIRONMENT&THRUSTERS&SENSORS) 

12.05 

215.03 

1.23 

8.39 


FLEXIBLE BODY 6 0.5 rpm ROLL RATE 


■ 

InWl 



ENVIR. (AERO+GRAV GRAD+CRYO) 


5.76 

1.04 

0.37 

THRUSTER ERRORS 

9.61 

9.54 

0.016 

7.39 

ROLL STAR TRACKER 
NOISE/QUANTIZATION 

0.29 

4.99 

0.55 

1.31 

SCIENCE TELESCOPE 
NOISE/QUANTIZATION 

8.08 

25.91 

0.14E-4 

0.35 

CONTROL GYRO DRIFT 

0.0050 

115.42 

0.47E-3 

0.0032 

CONTROL GYRO NOISE 

1.37 

17.62 

0.0015 

o 

b 

CD 

a 

CONTROL GYRO QUANTIZATION 


6.39 

0.002 

0.71 

DRAG FREE SENSOR NOISE 

0.029 

5.01 

0.35E-4 

3.18 

DRAG FREE SENSOR QUANTIZATION 

0.002 

5.03 

0.36E-5 

1.26 


The results summarized in this table indicate that the flexible body effects during 
periods where the spacecraft is stabilized at a constant roll rate (non transient periods) 
produce only small differences from the rigid body results. This indicates that the 
flexible body effects are of minimal importance to spacecraft control during quiescent 
periods. 
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Figure 3-1 LOS Error (.3 rpm, GSV and GSI, Environmental Dist Only) 


GV LOS ERROR (0.3 rpm), ENVIRONMENTAL DISTURBANCES ONLY 



Gl LOS ERROR (0.3 rpm), ENVIRONMENTAL DISTURBANCES ONLY 
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Figure 3-2 Roll and Translational Error (.3 rpm. Environmental Dist Only) 


ROLL ERROR (0.3 rpm), ENVIRONMENTAL DISTURBANCES ONLY 



TRANSLATION ERROR (0.3 rpm), ENVIRONMENTAL DISTURBANCES ONLY 
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Figure 3-9 LOS Error (.3 rpm, GSV and GSI, Control Gyro Drift) 
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Figure 3-10 Roll and Translational Error (.3 rpm, Control Gyro Drift) 


TRANSLATION ERROR (0.3 rpm), CONTROL GYRO DRIFT 
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Figure 3-12 Roll and Translational Error (.3 rpm, Control Gyro Noise) 
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Figure 3-13 LOS Error (.3 rpm, GSV and GSI, Control Gyro Quantization) 
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Figure 3-17 LOS Error (.3 rpm, GSV and GSI, Drag Free Sensor Quantization) 
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Figure 3-18 Roll and Translational Error (.3 rpm, Drag Free Sensor Quantization) 



TRANSLATION ERROR (0.3 rpm), DRAG FREE SENSOR QUANTIZATION 
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Figure 3-20 Roll and Translational Error (.3 rpm, All Disturbances and Errors) 
ROLL ERROR (0.3 rpm), ALL DISTURBANCES/ERRORS 
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Figure 3-21 LOS Error (.5 rpm, GSV and GSI, Environmental Dist Only) 


GV LOS ERROR (0.5 rpm), ENVIRONMENTAL DISTURBANCE ONLY 
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Figure 3-22 Roll and Translational Error (.5 rpm, Environmental Dist Only) 


ROLL ERROR (0.5 rpm), ENVIRONMENTAL DISTURBANCE ONLY 



TRANSLATION ERROR (0.5 rpm), ENVIRONMENTAL DISTURBANCE ONLY 
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Figure 3-25 LOS Error (.5 rpm, GSV and GSI, Roll Star Tracker Noise/Quantization) 
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Figure 3-33 LOS Error (.5 rpm, GSV and GSI, Control Gyro Quantization) 
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Figure 3-34 Roll and Translational Error (.5 rpm. Control Gyro Quantization) 
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Figure 3-36 Roll and Translational Error (.5 rpffi, Drag Free Sensor Noise) 
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Figure 3-37 LOS Error (.5 rpm, GSV and GSI, Drag Free Sensor Quantization) 
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3.1. Thrust Requirements 

The total thrust / flow requirements were examined for each of the sensor and 
actuator errors to determine which had a significant effect. In each case, the 
aerodynamic, gravity gradient and magnetic disturbances were applied. The 2000 
second run in each case started during the Guide Star Invalid period and switched to 
Guide Star Valid phase at approximately 900 seconds. A plot was produced for each 
case which included total thrust required, control thrust required, and flow (null dump) 
thrust required. The flow command was a constant 6.5 mg/sec (8.3 mn thrust) with an 
l,p of 130 seconds. Figure 3.1-1 shows the thrust command with perfect sensors and 
actuators. The control thrust command varies between 4 and 7 mn, below the 8.3 mn 
flow command. Figure 3.1-2 shows the thrust command with a roll star tracker noise 
whose standard deviation is 5 arcsec. The control thrust command exceeds the 8.3 mn 
flow command in this case. 

Figure 3.1-3 shows the thrust command with a science telescope noise whose 
standard deviation is 22.36 marcsec. The control thrust command is very nearly the 
same as the perfect case. Figure 3.1-4 shows the thrust command with a control gyro 
noise whose standard deviation is 0.002236 arcsec/sec. The control thrust command 
is also nearly the same as the perfect case. Figure 3.1-5 shows the thrust command 
with a drag free sensor noise whose standard deviation is 3.354 nm. The control thrust 
command varies between 4 and 8.4 mn, rarely exceeding the 8.3 mn flow command. 
Figure 3.1-6 shows the thrust command with a thruster noise whose standard deviation 
is 0.0559 mn. The control thrust command varies between 4 and 8 mn. Figure 3.1-7 
shows the thrust command with a roll star tracker quantization whose least significant 
bit weight is 0.5 arcsec. The control thrust command is again approximately the same 
as the perfect case. Figure 3.1-8 shows the thrust command with a science telescope 
quantization whose least significant bit weight is 0.25 marcsec. The control thrust 
command is again nearly the same as the perfect case. Figure 3.1-9 shows the thrust 
command with a pitch/yaw control gyro quantization whose least significant bit weight is 
1.3 marcsec. Once again the control thrust command is very nearly the same as the 
perfect case. Figure 3.1-10 shows the thrust command with a high range roll control 
gyro quantization whose least significant bit weight is 375 marcsec. The control thrust 
command exceeds the 8.3 mn flow command in this case. Figure 3.1-11 shows the 
thrust command with a low range roll control gyro quantization whose least significant 
bit weight is 37.5 marcsec. The control thrust command varies between 4 and 8.4 mn, 
rarely exceeding the 8.3 mn flow command. Figure 3.1-12 shows the thrust command 
with a drag free sensor quantization whose least significant bit weight is 1 .0 nm. The 
control thrust command is again nearly the same as the perfect case. Figure 3.1-13 
shows the thrust command with a thruster quantization whose least significant bit 
weight is 0.0025 mn. Once again little deviation from the perfect case is observed. 
Figures 3.1-14 through 3.1-16 show the thrust command with thruster hysteresis of 
0.100, 0.010, and 0.005 mn, respectively. The control thrust command in Figure 3.1-14 
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for 0.100 mn thruster hysteresis exceeds the 8.3 mn flow command, while the lower two 
hysteresis cases in Figures 3.1-15 and 3.1-16 are similar to the perfect case. Figure 
3.1-17 shows the thrust command with a control gyro drift of 0.003 arcsec/sec. The 
control thrust command again deviates little from the perfect case. 

These results indicate the potential sensor errors that can result in excessive 
thruster flow. These are the roll star tracker noise, roll control gyro quantization, 
particularly in the high range, and large thruster hysteresis, such as 0.1 0 mn. 
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Figure 3.1-9 Thrust Command For Pitch/Yaw Control Gyro Quantization 


THRUST COMMAND; PITCH/YAW CG QUANTIZATION 



THRUST COMMAND; ROLL CG QUANTIZATION 
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QPB Error Budget Evaluation Results 

The GPB error budget evaluation results from the GPB TREETOPS simulation are 
broken down into the Lockheed specified budget components. Table 4-1 shows the 
fine pointing, pitch/yaw, Guide Star Valid error budget components. Table 4-2 shows 
the fine pointing, pitch/yaw, Guide Star Invalid error budget components. Table 4-3 
shows the roll error budget components. Table 4-4 shows the translation error budget 
components. The short and long term errors as shown in the Lockheed error budget 
document* are combined in these tables. The errors not addressed in this study are left 
blank. The simulation errors quoted in the DETERMINATION, CONTROL, and 
summaiy blocks are taken from a simulation run where all of the subordinate errors are 
active. The number in parenthesis in each of these blocks is the RSS combination of 
the errors noted in the lower level blocks. Notice that the RSS totals and the composite 
simulation runs show approximately the same errors, indicating that the errors are 
essentially uncorrelated in producing control error. The Table 4-5 compares the 
observed control errors with the top level GPB control requirements and with the 
Lockheed budgeted allocations. Notice that the control errors predicted by the 
TREETOPS simulation are less than both the requirements and the Lockheed 
budgeted values. 
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Table 4-1 Fine Pointing Guide Star Valid Error Budget 

GS VALID POINTING ACCURACY 
Period < Orbit 

12.05 macs (RSS = 12.54 macs) 
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Table 4-2 Fine Pointing Guide Star Invalid Error Budget 

GS INVALID POINTING ACCURACY 
Period < Or bit 

215.0 macs (RSS = 197.3 macs) 
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Table 4-3 ROLL Error Budget 

~ ROLL ACCURACY 
Period < Orbit 

1.23 acs (RSS = 1.10 acs) 
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Table 4-4 Translation Error Budget 
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Table 4-5 GPB Simulation / Control Requirements Comparison 



Requirement 

Lockheed Budget 

Maximum Observed 
Error 

LOS Error - GSV 

20 mas (1 o) 

19.9 mas (1 o) 

12.05 mas (1 o) 

LOS Error - GSI 

2000 mas (1 o) 

56.9 mas (1 o) 

215 mas (1 a) 

Roll Error 

10 arcsec (1 o) 

1 0 arcsec (1 o) 

1 .23 arcsec (1 o) 

Translational 

Error 

100 nm (1 o) 

21 nm (1 o) 

8.39 nm (1 o) 
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5. Conclusions 

Models of fluid slosh and spacecraft flexibility were developed and added to the 
standard rigid body GPB simulation using the TREETOPS simulation package. 
Detailed slosh data was used to model the liquid center of mass motion produced by 
the liquid helium. A detailed flexible body model of the GPB spacecraft was obtained 
from Lockheed, processed, and integrated into the TREETOPS simulation structure. 

Performance results from the TREETOPS GPB simulation indicate that during 
quiescent periods of spacecraft activity, flexible body and slosh effects add only a small 
amount of additional control error over the standard rigid body response. The sloshing 
elements and flexible spacecraft components do not respond significantly to the 
spacecraft disturbances. It was discovered, however, that thruster flow limits are 
sensitive to specific sensor errors, such as roll star tracker noise, roll control gyro 
quantization, particularly in the high range, and large thruster hysteresis. If these errors 
are not filtered or held below specific values overcommanding of the thrusters is 
possible. 

A control error budget analysis was conducted and showed that all control 
requirements of the GPB program are met by the current configuration with significant 
margin. 
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